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Abstract 

Geometric phase that manifests itself in number of optic and nuclear experi- 
ments is shown to be a useful tool for realization of quantum computations in so 
called holonomic quantum computer model (HQCM) . This model is considered 
as an externally driven quantum system with adiabatic evolution law and finite 
number of the energy levels. The corresponding evolution operators represent 
quantum gates of HQCM. The explicit expression for the gates is derived both 
for one-qubit and for multi-qubit quantum gates as Abelian and non-Abelian 
geometric phases provided the energy levels to be time-independent or in other 
words for rotational adiabatic evolution of the system. Application of non- 
adiabatic geometric-like phases in quantum computations is also discussed for 
a Caldeira-Legett-type model (one-qubit gates) and for the spin 3/2 quadrupole 
NMR model (two-qubit gates). Generic quantum gates for these two models 
are derived. The possibility of construction of the universal quantum gates in 
both cases is shown. 

Keywords: Quantum computer, Berry phase, Non-adiabatic geometric phase, 
Two-qubit gates 

1 Introduction 

The conceptions of quantum computer (QC) and quantum computation developed 
in 80-th [T], [2] were found to be fruitful both for computer science and mathemat- 
ics as well as for physics [3]. Although a device being able to perform quantum 
computations is now far away from practical realization, there is a great number 
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of theoretical proposals of such a construct (see e.g. [6]-[T8]). Intensive investi- 
gations on quantum information theory (see e.g. [3], [5] for a reference source on 
this subject) refreshed some interest on Berry phase [56]. The idea of using unitary 
transformations produced by Berry phase as quantum computations is proposed in 
|19j . |20j and first realized in [21], [U] in a concrete model of holonomic quantum 
computer where the degenerate states of laser beams in non-linear Kerr cell are 
interpreted as qubits. For other references where Abelian Berry phase is considered 
in the context of quantum computer see e.g. |22j - ^25j. If the corresponding energy 
level is degenerate non- Abelian phase takes place |57j that is actually a matrix mix- 
ing the states with the same energy. For further references on quantum computation 
based on non- Abelian geometric phase see e.g. [47]-|53j. 

On the other hand non-adiabatic analogue of Berry phase can exist and be mea- 
sured if transitions in a given statistical ensemble do not lead to loose of coherence 
|71j . For loose of coherence in quantum computations related to geometric phase 
see [1I]-[16]. Thus it is also possible to use the corresponding unitary operators to 
realize quantum gates. This fact has been noticed in [26], [27]. After that a lot of 
papers was published where the non-adiabatic phase is applied to realize the basic 
gates in different models of QC such as different NMR schemes |28j-|35j. ion traps 
[36] , [37] . quantum dots [38], [39], and superconducting nanocirquits [30] . 

To analyze a concrete scheme for quantum computation based on geometric 
phase it is desirable to be aware of analytical expression for the evolution operator of 
the system at least at the moment when the measurement is performed. This article 
is concentrated on the computational aspect of geometric phase for the models which 
are relevant to QC. It should be emphasized that the form of the expression for the 
phase and the possibility of the derivation of such a formula itself thoroughly depend 
on the group-theoretic structure of the corresponding Hamiltonian. Therefore a 
method of the geometric phase calculation which would be more or less universal at 
least in the adiabatic case can appear to be useful. The material is divided in two 
parts. In section [2] the adiabatic geometric phase is considered. In subsection 1 2 . 1 1 we 
analyze the difficulties appearing in calculation of the Abelian adiabatic geometric 
phase (Berry's phase) and propose a method of its explicit derivation for the case 
of the symmetric time-dependent Hamiltonian with constant non-degenerate energy 
levels. The symmetry of the Hamiltonian is supposed to reduce the Hamiltonian 
to that of a system with finite number of energy levels. In subsection 12.21 this 
method is generalized for the case when degeneration is present. In section [3] we 
consider non-adiabatic phase which can only conditionally be called "geometric" for 
its dependence on concrete details of the dynamics. For this reason it is not possible 
to work out more or less general approach to the calculation of the non-adiabatic 
phase. Therefore two concrete cases are considered. In subsection 13. II application of 
the Abelian non-adiabatic phase to one-qubit computation in a Caldeira-Legett-type 
model is cosidered. In subsection l3.2l we present an example of both non- Abelian and 
non-adiabatic phase computation in spin-3/2 quadrupole NMR resonance model. 
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2 Adiabatic Geometric Phase 
2.1 Abelian Berry's Phase 

Here we consider a possible method of the adiabatic phase computation that seems 
to be effective in a broad range of practically relevant cases. Berry phase is a conse- 
quence of the adiabatic (or Born-Fock) theorem [55] which states that a parametric 
quantum system depending on a set of slowly (adiabatically) evolving parameters 
Ri{t), i = 1, . . . N behaves in a quasi-stationary manner 

H{R)\n{R)>= En{R)\n{R)>, R = {Ri, . . . , Rn) (1) 

where H{R) is the corresponding Hamiltonian and no energy level degeneration is 
assumed. The adiabaticity condition means that the frequencies = En{R)/h 

are much greater than the characteristic Fourier frequencies of Ri{t). Thus the 
eigenvectors \n{R) > evolve like 

\n{R) >= S{R)\no >, |no >= HR{0)) >, SS^ = 1 (2) 

with unitary rotation S describing the natural variation of \n(R) > due to that of 
R(t). It corresponds to the following evolution law of the Hamiltonian 

H{t) = S{R)Hoit)S\R) (3) 

where Ho{t) is diagonal in the basis {|no >}. What is the solution of the non- 
stationary Schrodinger equation 

in^^\m>=H{t)\m> (4) 

for this case? A natural hypothesis would be that the evolution operator for 1-0 > 

\m >= u{t)\m > 

has the form 

U{t) = S{R)^{t) 

where S is defined by ([2]) and ^> simply produces the dynamic phase 



6(^)|n(i^(^) >= exp j -i/h j En{T)dT j \n{R{t) > . 



(5) 



It is based on the analogue with the stationary case where evolution is simply rep- 
resented by the dynamical phase factor exp(—i/hEnt). Berry first observed |^56] 
that the hypothesis is wrong. To see this it is sufficient to represent U in the form 
U{t) = S{R)V{t) and substitute it into the Schrodinger equation 

ih^U{t)=H{t)U{t). (6) 
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It gives 

ih{VV^ + S^VrSR) = Ho{t). (7) 

Now one can see that V cannot be simply $ because it has to cancel the second 
term in the right hand side of d?]) besides of -f^o- It follows from ([TDthat 

v{t) = fit)m 

where is determined by ^ and the following equation is valid for T{t): 

f f t |no >= -{S^VrS)R \no > . (8) 

It results in the evolution law for the state vector corresponding to the n-th energy 
level 

>= e-V'^^-W e^^"W \n{R{t)) > (9) 
where ^n{t) is the phase factor in the right-hand side of ([5]) and 7„(f) is given by 
t 

^n{t) = j An{R{T))R{T) dr, A(l?) = i < n{R)\ VRuiR) > . (10) 



Phase 7n(t) becomes purely geometric while R{t) evolves cyclically: R{T) = R{0) 

-fr,{T) = 7„(C) = dR. (11) 

c 

Here the integration contour C is a closed curve in the parameter space described 
by R{t) as a radius- vector. It is easily seen from (fTT]l that 7„(C) does not depend 
on the concrete details of the system's dynamic if the adiabatic condition is held. 

In this article we are interested in computing of 7n(C) in the most general case. 
The problem of derivation of 7ri(C) was solved in various particular cases in large 
number of articles some years ago. First we would like to note that the straightfor- 
ward formula 

R ^ ^ . V- {yRH{R))mn X {VRH{R))nm 

= {E^{R) - {Em{R)Y ^''^ 

derived by Berry [56] by making use of the identity 

< mlV/^ln > 
< m|Vn >= — — — — , m ^ n 

[En — Em) 

has not (despite of it's beauty) much practical use because to apply it one should 
establish the analytical dependence of all En on R that is not a realistic task exclud- 
ing some special cases. To see this one should attempt to apply formula (fT2l) to the 
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case of a 3- level system substituting a generic solution En{R) of the corresponding 
cubic equation therein. 

It was first noticed in ^Ij that symmetries of the Hamiltonian H(R) play an 
important role in computing of 7^- Indeed if one represents the result of the periodic 
motion H{0) = H{T) as 

|V(r) >= U{T)\m > (13) 

where as it follows from ([1]) U{T) must commute with H{0) so it can be represented 
as an exponent containing a linear combination of operators which must com- 
mute with H{0) as well. Thus the operators are integrals of motion and describe 
certain symmetries of the given system. Therefore in what follows we restrict our- 
selves with such systems whose Hamiltonian is an element of a finite Lie algebra. 
This assumption immediately gives the group-theoretic structure of U{T): 



UiT) = exp{iya,Hi] (14) 




where Hi are all linearly independent elements of the Cartan subalgebra and Oj are 
some coefficients. Thus the problem reduces to computing of the coefficients Oj. In 
the simplest case of Lie algebras consisting of three elements this problem can be 
easily solved [62], [6l|, [65] for physically relevant cases of Heisenberg-Weyl algebra, 
su{2) and sm(I.I). In each of them the evolution operator has the form 

U{t) = exp (c(t)X+ - C{t)X^) exp 0(^X3) (15) 

provided the initial Hamiltonian is proportional to X3 where X± and X3 are the 
corresponding generators of the algebras above. Their expressions for each concrete 
case are given in table [U The Hamiltonian for su{2) case describes an arbitrary 
spin in the magnetic field so all J's are the angular momentum operators: J± = 
l/2(Ji lb J2). sn(l.l) case corresponds to the evolution of squeezed states [66] of 
light in non-linear optics. Here = /2, K- = c? 12 and Ky, = a"'"a + 1/2 where 
a, are usual bosonic annihilation and creation operators. The last case represents 
a harmonic oscillator interacting with the time-dependent electric field. The simple 
commutation relations in these three algebras admit direct computation of Berry's 
phase [S2], [Ml, [65]. 

-fm = m ^uj{0 = j dAioiO (16) 
c s 

where m is an eigenvalue of the corresponding X3, S is the surface in the parameter 
space bounded by the closed curve C and the expressions for u){^) and its exter- 
nal derivative d A uj{(^) are given for each case in table [21 The geometric sense of 
the derived phase factor is the integral curvature over the surface bounded by the 
contour C on the manifold the evolution operator belongs to. This manifold can 
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Table 1: Expressions for the operators X±, X3. 



Alg(;l)ra 


x+ 




X, 


Commutators 










[J3,j±] = ±i± 


su{2) 


J+ 




J3 


[J+,J_]=2J3 




















stt(l.l) 


K+ 


k_ 


K3 












[k+,k_] = -2k3 










[i,a+] = [i,a]=0 


H-W 


a+ 


a 


i 










[a, a+] = i 



be generally expressed in the form G/H where G is the group manifold and H is 

that of the stationary subgroup, i.e. the group whose Lie algebra consists of all 
operators commuting with H{0) (in these three cases it is always U{1)). It is sphere 
in the case of su{2), two-sheet hyperboloid in the case of su(l.l) and plane in the 
case of Heisenberg-Weyl group. To complete the computation one has to establish 
correspondence between the complex parameter ^ and physical parameters of the 
Hamiltonian. Let us do that for SU{2). It is worth to notice that the result is 
completely determined by the geometric properties of the group and does not de- 
pend on the concrete representation. For this reason one can chose the fundamental 
representation of SU{2) to simplify the derivation. Thus we take H{<S) = lob era 
which corresponds to the initial eigenvectors |± >= (1(0), 0(1))^ (T denotes trans- 
position). The evolution operator generally parametrized by the spherical as 



cos- 



sin| e-'^\ 



sin-| e'^ cos^ 



(17) 



rotates H{0) into H{t) = LOBntr where n = (cos (/3sini9/2, sin99sin'!9/2, cost?/2) 
determines the direction of the magnetic field. On the other hand the direct com- 
putation of exp((^J+ — C*J-) gives for this representation \(\ = i?/2, arg(^ = ip + 
It leads to the well known expressions for the fictitious "strength field" 



B 



± 



^ 2 i?3 



(18) 



where R denotes the true magnetic field vector in order not to confuse it with the 
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Table 2: Expressions for the forms uj and d Auj. 



Algebra 




d A w(^) 


Relation to Q 


su{2) 


CdC-CdC 


2d^AdC 


iei=tan(|C|), 
arg ^ = arg C, 


sn(l.l) 


CdC-CdC 


2dCAdC 


lei =tanh(|C|), 
arg ^ = arg C, 


H-W 


CdC - Cd^ 


2d^AdC 


e = c 



fictitious one which determines the resulting Berry's phase. It should be noted that 
the correspondence i? — > realizes the stereographic projection of the sphere with 
the coordinates i?, (/? on the plane that points are labeled by ^. The other two cases 
of SU{1.1) and Heisenberg-Weyl groups can be considered in a similar manner. 

Re-derivation of these simplest results has the intention to extract a universal 
idea of computing the geometric phase in more or less general case. For the sake of 
certainty let us assume the symmetry algebra of the Hamiltonian to be semisimple. 
It means that the generic evolution operator can be represented in the form 

u{t)= n (19) 

where A_|_ denotes the set of the positive roots a and each Ua is analogous to (fT5]) 
(see also tabled for su{2) and s'u(l.l) cases): 

?7„(t) = exp (a(t) - C{t) E-a) (20) 

where the standard notations for the Cartan basis [67J 

[Hf3,E±a\ = ^a{Hp) E±a [Ea^E-a] = Ha, [Ea, Ef^] = N^pEa+p (21) 

are used. The pairs of generators E±a are analogous for J± in su{2) and i^^'s are 
that of J3. Here a{Hp) and N^fi are constants that can be chosen rational and 
integer correspondingly. Taking account of the consideration above leads to some 
more detailed form for one-cycle evolution operator U{T) (jl4p 

IJ{C) = exp j i a''{C)Ha ) . (22) 

\ "GA+ / 
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Each pair {Ea,E^a) makes besides of the trivial group-theoretic contribution Ha 
which produces the corresponding quantum number also a non-trivial one reflecting 
adiabatic dynamics of the system 



a"(C) = 




c 



where 9a generally depends on all Ca{t)- Thus to solve the problem one has to find 
this dependence making use of commutation relations (j2ip and then establish the 
connection between the parameters Ca(^) ^^id the natural set of parameters R of 
the Hamiltonian. Unfortunately the hope to obtain a solution of even one of these 
two tasks that would be a non-trivial generalization of the above examples is not 
realistic. Neither the first part of the problem nor the second one could be solved in 
a way resulting in physically relevant explicit formulas having practical use. First 
the 1-forms 6a fulfill Maurer-Cartan equations that express the quantity W dU in 
terms of the 1-forms oja and Oi 

U\C) dU{C) = iiuJaiC) Ea + d,{C) Hi) (24) 

where the index i labels all linearly independent generators of the Cartan subalgebra 
(not all Ha are so). For commutation relations (j2ip these equations take the form 

d A = A + 1/2 ujf^ A uj^ (25) 

d^e'' = 1/2 Ciy^ ujf^ A uj^ (26) 

Equations (|25p . (I26p describe the parallel transport on the coset manifold G/H. 
The possibility to solve them depends on the manifold's symmetry and of course 
is entirely determined by the structure constants C.. that are built from the root 
vectors a{Hp) and the constants Nap- The general solution of this system can 
be constructed for very high symmetry of symmetric spaces |68] where the whole 
algebra can be split in two subsets X and Y such that 

[Y,Y\CY, [Y,X]CX, [X,X]CY. 

It is seen from (I2ip that the last condition is generally speaking not valid for our 
case because not all A''^^^ are zeroes. Its geometric sense is that the considered coset 
spaces G/H are of more general symmetry type than symmetric spaces. Thus for 
G = SU (n) the space 

SU{n)/U{l) X U{1) ... X U{1) (27) 

^ V ' 

n—1 times 

belongs to the more general class of Kahlerian spaces. The general solution of (j25p . 
(j26p for the types of spaces we are interested in is not obtained so far. Therefore 
the practical use of these equations is not high. Moreover the solution of the second 
part of the problem discussed is not possible for the same reason. 
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A simple and effective metfiod of practical computation of geometric phase where 
it is not necessary to find the forms w", 9^ is proposed in [75j. For this purpose we 
have to make some assumptions. First we regard the Hamiltonian to belong to 
a finite irreducible representation of a semisimple Lie algebra therefore Hq in ([3|) 
can always be represented as a finite matrix Hq = Ri,{t) Hi where the set {Hi} 
is a basis of the Cartan subalgebra and Ri{t) are parameters. Then we suppose 
the energy levels Em to be constants. It corresponds to a rotation-type evolution 
([3]) where Hq does not depend on t. Such a situation takes place practically in all 
experiments on the geometric phase measurement. This makes it possible to regard 
Em as additional secondary parameters to be found just once (may be numerically). 
The third assumption is that the spectrum remains always non-degenerate i.e., no 
crossing of energy levels occurs. As the spectrum of the Hamiltonian is finite, the 
state vector \ipm) is a unit vector m in C", so Am is 

Am = - {m*dm — mdm*). (28) 

As the evolution is adiabatic, the spectrum of H{t) remains always non-degenerate 
if it was so at the initial time. Then there is always a nonzero main minor of H — Em 
which we assume to consist always of the first n — 1 lines and columns oi H — Em ■ 
Denoting the matrix consisting of the first n — 1 lines and columns of H by H±^ we 
come to the condition 

det(i/x - ^m) / (29) 
Making use of this condition one can represent n in the uniform coordinates 

and express $,rn in terms of Hij for 1 <i,j <n — 1 and Em- 

= (-f^-L ~ Em)~^ h, hi = —Hin, (30) 

where ^ is a vector in C""^ but not in C". Thus we have for Am 

where ^m is completely determined by ([30]) . Note that the result obtained is purely 
geometrical because it can be expressed of the Kahlerian potential 

F = log(l + |z|2) 

where z is a vector in consisting of n(n — l)/2 independent components of all 
^m- The function F{z,z*) determines all the geometrical properties of the state 
space (|27p . Particularly its metric tensor is 

_ d'^F{z,z*) 
~ dz,dz* ■ 
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It should be emphasized that the simpHfication of the problem reached here is based 
on the fact that Em = so one can include it in new parameters and use them rather 
than Hi. Therefore the dependence of Em on R is not required. One can calculate 
Em numerically and substitute it into the formulas regarding this quantity as one 
more external parameter. Moreover to find 7^ one needs only the energy Em but 
not the whole spectrum as in (jl2p . It can become an important issue if one considers 
partially solvable models. The requirement Em = is sufficient because otherwise 
one has to solve the secular equation at each moment t that is equivalent to the 
numerical solution of the non-stationary Schrodinger equation itself and therefore it 
makes the discussed method useless. 

Let us now consider some simple applications of the proposed method. First 
let us see how it works for the trivial case n = 2. ([1]) reduces then to two linearly 
dependent equations 

{B3TB)^+iBi-iB2) = 
{B,+iB2)C + {-B3TB) = 

Here we returned to the usual notations of the magnetic field components Bi and 
zizB = ±|S is the energy of the state |ib >. Choosing one of the equations and 
taking the spherical coordinates we come to one of the relations 

^ = -tani?/2 e"*'^ ^ = cot ^9/2 e"^'^. 

for the upper and lower sign correspondingly. Thus these are the coordinates of 
stereographic projection made from the north (south) pole of the sphere. Substitut- 
ing it into the formula for uj±{(,) (see Table [2]) and integrating over a contour C we 
get the well known result [56j 

^^ = ^2} i + |gp = ^2 

c 

where ri(C) is the solid angle corresponding to the closed contour C on the sphere. 

One more example which is less trivial is a generic three-level system. The 
k-th eigenvector is then two-dimensional and some trivial algebra gives for its 
components 

a = Ai/Ao, 6 = A2/Ao, 

Aq = {Hii — Ek)[H22 — Ek) — |-ffi2p 

(32) 

Ai = H23H12 — His{H22 — Ek) 

A2 = — H2'i{Hii — Ek) 

Here we have omitted where possible the index k. Substitution of these expressions 
into (j3ip gives the final formula for this case. Note that the use of formula (jl2p 
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here would lead to sufficient computational difficulties even after making further 
simplifying assumptions [69j. The proposed approach makes concrete calculations 
visibly easier and more compact although the final formulas are not of esthetic value. 
For the illustrative purpose we take here the case when all Hij, i < j but H12 do 
not depend on t. Then substitution of (f32]) into (fSTI) gives 

Wfc(C) = iCk [A- Dk sm{(pi2 + (p23 - (Pn)] d(j)i2, (33) 



(-^ ^ 1-^^1311^2311^121 

*^ Al{Ek) + \A,{Ek)\^ + \A2{EkW' ^ ^ 



A = l/|//i3|2-l/|i/23|', (35) 



Dk = Hn + H22-2Ek (36) 

where (pij are arguments of the complex numbers Hij. As it was discussed above 
the condition Ao(t) 7^ is supposed to be held everywhere on C. For other cases of 
the geometric phase in the 3- level system see |63j . 

2.2 Non-Abelian Wilczek— Zee phase 

Now we proceed with a more general case of degenerate spectrum. Quantum com- 
putation for this case generated by an adiabatic loop in the control manifold is de- 
termined by the same quasi-stationary Schrodinger equation ([1]) where each energy 
level Em corresponds to a set of eigenstates \ma>, a = 1, ... dm- Cyclic evolution 
of the parameters results in 

\ma{T)>= Uab{T) |mb(0)> (37) 
where the matrix?/ is presented by a P-ordered exponent 



U{e) =Pexp A, 



{Am)ab = i <'mb\ dma > . (38) 



In this section we generalize the proposed approach to the geometric phase compu- 
tation for the generic case of degenerate energy levels |75| . 

The set of eigenvectors $ma^ ^ — ^1 ■■■idm must obey the equation 

{Hf-^ -Em)ima = hCa. (39) 

Here the matrix H^f"^^ is constructed from the first n — dm lines and columns of H, 
Ca are arbitrary dm-dimensional vectors and h is the following {n — dm) x dm-matrix: 

Hi 



h 



\Hn-dm,n-dm+l ■ ■ ■ Hn-dm,n j 
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Of course it has sense only if the condition 



(40) 



is vahd along the evolution process. The set of vectors must be orthogonalized 
by the standard Gram algorithm and after that we get the orthonormal set of the 
eigenvectors Za (here and below we has omitted the index m) in the form 



1 



det Ta- 



Xi\ 



a-l 



(41) 



where = (^j,, q) and is chosen to be the standard orthogonal set 

a 

Ca = (O.-. 1 ...0). The matrices Ta are determined by 



1 + ZiZa, 



(42) 



where the (n — dm) x a-matrix Za consists of a first lines of the (n — dm) x dm-matrix 
Z = {H^'"^ — Em)~^ h. Using (j42p and (|4ip we come to the final expression for the 
matrix- valued 1-form A: 



2 det(l + zl_^Za-i) det(l + 4~i^fe-0 ' 



l<i<a, 1 < J < 6, 



(43) 



where 



9ab 



-pi p*J 
^ a'- b ^ 



U^ab = i^j C?Im(5ab) ^i) + 



min(a,b) 

E 



dlraig'Jb), 



and is the cofactor of in Tq. Note that the change of our basis by = 
Uab{^) Cb leads to a standard gauge transformation of A 

A' = UAU^ + i{dU)UK 

The formula (|43p is the desired expression of A in terms of the matrix elements 
of the Hamiltonian. It is correct if condition (j40p is valid. It is not nevertheless 
a principal restriction because dm does not depend on time due to adiabaticity of 
the evolution and there is always at least one nonzero n — dm-order minor of H. 
Then, if the minor we choose vanishes somewhere on the loop C one can always 
take local coordinates such that the techniques considered is applicable on each 
segment of C. It should also be noted at the end of this section that the idea of 
physical realization of the quantum gates based on the concrete system driven by 
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external electromagnetic fields appears if one takes into account that for An 
can be realized by means of ordinary bosonic creation and annihilation operators, 
namely = Oifoj for some 1 < i, j < n. Then Eq, represents nothing but two-mode 
squeezing operator. Thus the model considered can be applied to optical HQC with 
n laser beams (the case n=2 is considered in [21]) and the logical gates Ua are just 
two-qubit transformations realized by transformation of two laser beams. 

The method presented here enables one to build in principal any computation 
for HQC described by a Hamiltonian with a stationary spectrum in terms of ex- 
perimentally measured values exactly the matrix elements of the Hamiltonian. The 
method depends weakly on the dimension of the qubit space which other models 
based on various parameterizations of the system's evolution operator are very sen- 
sitive to. Application of this method to a concrete physical model will be discussed 
elsewhere. 

3 Non-Adiabatic Geometric Phase 
3.1 Abelian Non-Adiabatic Phase 

The adiabatic condition of quantum system's evolution is strong enough to restrict 
sufficiently the scope of the search for realistic candidates for practical realization 
of quantum computations despite of some attractive features of the adiabatic case 
such as fault tolerance due to independence of the evolution law on the details of the 
parameters' dynamics etc. Therefore it is desirable to find physically relevant cases 
for which on the one hand this condition would be not necessary but on the other 
hand the coherency in such a system would be not yet violated so that the notion 
of the phase shift itself could have physical sense. As the adiabatic theorem is no 
longer valid the property of universality of the system's dynamics (independence of 
the concrete form of the functions Ili{t)) is no more preserved and the evolution 
law is sufficiently more complicated. Then one cannot hope to carry out a general 
approach to derivation of the corresponding phase shift because in each case it 
depends on the fine details of the parameters variation. For the same reason the 
phase can be called "geometric" only conditionally because geometric intuition is 
no more helpful for this case e.g. the result can be represented as an integral over 
t rather than over a contour that expresses mathematically the thesis above. On 
the other hand non-adiabatic conditional geometric phase that was theoretically 
predicted in [60] can be measured if transitions taking place in the system do not 
lead to decoherence [TT]. Therefore it is also possible to use the corresponding 
unitary operators to perform quantum calculations. This fact has been noticed in 
[26], [27] (see also |28]-|35| for further references). 

Let us consider a parametric quantum system described by the Hamiltonian 
H{R), where R{t) is a set of arbitrarily evolving parameters. We suppose that 
evolution of the Hamiltonian is determined by unitary rotation ([3]) Looking for 
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particular solutions of the Schrodinger equation dH) (here we supposed ?i = 1) we 
take a rotating frame by assigning -(/^(t) = U{t) ip{t) and get in such a way 

i^-l = {H,{t)-^uHt)U{t))m■ (44) 

Of course, this transformation generally does not help to solve equation Q due to 
the fact that the algebraic structure of the coupling term —i WU can appear to be 
rather complicated and the last generally does not commute with Hq. However if 
a receipt is known how to evaluate the last term in ()44p . further solution of this 
equation is straightforward: 

^(t) = e-^<^" W Texp (^-i U\t)U{t) dr^ V'(O), (45) 

where (l)n{t) = Jq En{T) dr is so called dynamic phase, T denotes time-ordering and 
En are elements of Hq that is by definition diagonal. Of course if there is no way to 
find U\t)U{t), expression ([^5]) is useless. 

Let us illustrate it for the simplest case of spin 1/2 in the non-adiabatically 
rotating magnetic field [72] . Uniform rotation in the plane 'd = const is represented 

by 

H{t) = e^*'^«* ^3 e-'^ ^2 g-ii? J2 ^Ti^Rt J3 (4g) 

where the sign + (— ) corresponds to the left (right) polarization. Application of 
(j44p to (I46p gives for the Hamiltonian in the rotating frame 



H,{t) = e-''^'^ (nJs) e-''''^±u:Rh. (47) 

To diagonalize Hamiltonian (j47p one has to apply one more rotation to it 

H2{t) = V Hi{t)V\ V = e^^*^^ 

where the angle t?* does not coincide with due to the second non-adiabatic term. 
It should rather fulfill the condition 

.a.r = ff (48) 

cos V ± UJr/U 

The second term in the denominator of (|48p is the measure of non-adiabaticity of 
the motion. It is clear that the angle -d* replaces the usual azimuthal angle in the 
formula for the geometric-like phase: 

7^ = zp ma 2tt (1 - cos •&*) (49) 



where ma is the third spin projection. Formula (j49p is a natural generalization of 
the usual Berry's formula for the adiabatic case and coincides with it in the limit 
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iVpi/Vt — > 0. Note that the dependence of the result on ujn reflects the fact the phase 
is no longer truly geometric because lor characterizes the rotation velocity and thus 
the velocity of the motion along the contour in the parameter space. 

As an example of the application of the non-adiabatic formula above we propose 
a realization of quantum gates for a concrete 4-level quantum system driven by 
external magnetic field [7^. Let us consider a system of two qubits in a bosonic 
environment described by the Hamiltonian 

H = Hs + Hb + Hsb, (50) 
where Hs is the Hamiltonian of two coupled spins 

Hs = + i/if* = a,i l2 + l2 » CT,2 + ^ <Jz2, (51) 

where J is the coupling constant, Hb is the Hamiltonian of the bosonic environment 

Hb = Y, ^bkibth + 1/2), (52) 
k 

and Hsb is the Hamiltonian of the spin- enviroment interaction. 

Hsb = H^sB + (53) 

H^SB = Si^^Y^^^-kbt + glkh) a = 1,2. (54) 

k 

Here 

= l2, 5p) = l2®a,2, 

Uz is the third Pauli matrix, I2 is 2 x 2 unit matrix, b'^,bk are bosonic creation 
and annihilation operators and Qak are complex constants. We assume that the 
two spins under consideration are not identical so that wqi 7^ ^02- The Hamiltonian 
determined by (j50p - (j54p is a natural generalization of Caldeira-Legett Hamiltonian 
|74j for the case of two non-interacting spins. Let such a system be placed in the 
magnetic field affecting the spins but not the phonon modes. The only change to 
be made in the spin part ([5T]) is the substitution 

i^sCTz — > Ba, 

Three components of B represent a control set for the qubits under consideration. 
Evolution of B{t) generates evolution of the reduced density matrix ps{t) that 
describes the spin dynamics 

i^PM = HsPs{t), Ps{t) = U{t)p{0)U+it). (55) 
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Thus given curve in the control space corresponds to a quantum calculation in 
which each qubit is to be processed independently. To obtain such a calculation as 
a function of control parameters we first recall some common issues of spin dynamics. 
We consider the external magnetic field as a superposition of a constant component 
and a circular polarized wave: 

B = Bo + -Bie^'^«*, (56) 

where Bq is perpendicular to -Bi. It is well known that the case of the circular 
polarization is exactly solvable. The evolution of an individual spin corresponding 
to the Hamiltonian 

H = -fiB (57) 

is determined by (fT5]) where X±, X3 are replaced by S± = S^^iSy, Sz correspond- 
ingly and 

({t) = \C{t)\e^p{iAu;t + ia{t) + iTr/2) , (58) 
uj± sin(J7t/2) 



m\ 



a{t) = arctan ^-^^ tan(r2t/2)^ , 

0(t) = _,^^(^^n2+6ni), (59) 

where Au = W|| — lur, = (Au;)^ + ^± and are Rabi frequencies corre- 
sponding to Bq and Bi respectively and finally n is the unit vector along Bi. 

It is known [71j that the pure states acquire within the rotating wave approxi- 
mation a phase factor that after one complete cycle T = 2tt/ujr is: 

\m{T) >= exp(— i0D + 17) |m(0) >, (60) 

where m is the azimuthal quantum number and the phase is split in two parts: 
dynamic 

= 27rm — cos{6 — 6*) 



and geometrical 



7 = — 27rmcos6'*, (61) 



where cos9 = Bq/B {B = Bq + Bi) and 9* is determined by formula (|49|) . The 
phase shift between the states | it 1/2 > results then in 



A(l)g = -2-Kcose* 



(62) 
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that is nothing but the sohd angle enclosed by the closed curve B{0) = B[T) on the 
Bloch sphere. If the rotation is slow such that uj^/il. then 6* ^ 9 and phase 
shift (|62p coincides with the usual Berry phase. 

Thus the adiabaticity condition is not really necessary for obtaining of the ge- 
ometrical phase in an ensemble of spins if the decoherence time is much greater 
than T. Therefore one can attempt to use this phase to get quantum gates such 
as CNOT. Calculation of the corresponding phase factors is rather straightforward 
because the free and the coupling parts of the spin Hamiltonian commute with each 
other 

'HP,Hf] =0. 

Therefore the coupling part can be diagonalized simultaneously with the free part 
by applying of the transformation U = Ui <S> U2 where are the diagonalizing 
matrices for each single-spin Hamiltonian respectively. This simple fact together 
with the following obvious identity 

U^U = ulUi O I2 + I2 ulu2 

the final formula for the part of the evolution operator that stands for the non- 
adiabatic geometric phase 

Ug = exp(— 27r« cos 9\ Siz) tX" exp(— 27rz cos 62 S2z), (63) 

where 

, n* sin 6*1 sin 6^2 
tanfc'^ = — — , tanfc'a = jrr- 

COS t)l + U>Ji/ill COSU2 + UJr/U2 

and 



cos 01 = Woi/^^l, ^1 = UJqi + to 



1) 



COs6'2 = ^^02/^^2, f^2 = ^^02 + '^1- 

Note that gate (j63|) is symmetric with respect to the spin transposition as it should 
be and does not depend on J that is typical for geometrical phase in spin systems 
where the phase depends only on the position drawn by the vector B on the Bloch 
sphere. As J does not affect this position, it is absent in the final result. We do not 
consider here the dynamic phase determining by the factor 

Ud = exj> [-\HsT 

It is so because one can eliminate it by making use of the net effect of the compound 
transformation proposed in [22]. After this transformation that is generated by two 
different specifically chosen contours the dynamic phase acquired by the different 
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spin states becomes the same and the geometric phase of each state is counted twice. 
After that we get (up to a global phase) the following quantum gate 



Ua 



gii-'yi+ii) 

V 



\ 



(64) 



-4(71+72) J 



Thus we have constructed the quantum gate, which possess the advantage to be 
fault tolerant with respect to some kinds of errors such as the error of the amphtude 
control of B. On the other hand this approach makes it possible to get rid of the 
adiabaticity condition that strongly restricts the applicability of the gate. Instead of 
this condition one needs some more weak one: r ^ ^r^^ where r is the decoherence 
time. 



3.2 Non-Abelian and Non-adiabatic Phase 

In this section we give an example of both non-Abelian and non-adiabatic phase for 
a concrete 4- level quantum system driven by external magnetic field [76] . Let us con- 
sider a spin-3/2 system with quadrupole interaction. Physically it can be thought of 
as a single spin-3/2 nucleus. A coherent ensemble of such nuclei manifest geometric 
phase when placed in rotating magnetic field. This phase is non-Abelian due to de- 
generate energy levels with respect to the sign of the spin projection. Depending on 
the experiment setup the phase can be both adiabatic as in Rb experiment by Tycko 
[70j and non-adiabatic as in Xe experiment by Appelt et al [71]. This non- Abelian 
phase results in mixing of ±1/2 states in one subspace and ±3/2 in another one 
and thus can be regarded as a 2-qubit gate. The gate is generated by a non-Abelian 
effective gauge potential A that is the subject of computation in this section. 

We assume the condition of the ^^^Xe NMR experiment to be held so one does 
not need to trouble about the coherency in the system. The last is described by the 
following Hamiltonian in the frame where the magnetic field is parallel to the z-axis 
{h=l) 

Ho = - + !))• (65) 

Here and in what follows we omitted the hat symbol over all J's for the sake of 
simplicity. For a spin-3/2 system we choose the third projection of the angular 
momentum in the form 



/3/2 








^ 





-3/2 














1/2 





V 








-1/2/ 



/3/2a3 

V 1/2^73 



(66) 
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Then two other projection operators are 



/ 



V3/2 

V ^3/2 1 



V3/2 \ 
V3/2 
1 

/ 







J2 



-iters' 



(67) 



(68) 



/ -\/3/2\ 

^/3/2 

V3/2 -i " ^^ as CT2 y 

V-\/3/2 i / 

In the laboratory frame the Hamiltonian takes the form 

H = uJo{{Jnf - l/3i(j + 1)) = e-'^^-'e-'^^^ Ho e''P^^e''^^\ 

Rotation around the z-axis means that (/? = ujit and one should perform the unitary 
transformation 



1^ >= Ui 1-0 >, 
In the rotating frame we get 



(wo-'a ~ '^i-'sje -7-, 



(70) 
(71) 



where 



Js 



jej2 



J?, e 



Expression (jTTj) is equivalent to 



' ujQ — -^uJi cos a (T3 



V3/2 



(72) 



uj\\f?>l2 — Wo — ^t^i cos (T3 + (Ji sin (Ti^ 

It is convenient to diagonalize this matrix in two steps. First we get rid of a\ in the 
last matrix element by performing of the block-diagonal transformation 



[/2 = diag(l,e-^"'^^), 
where tana = 2tan0. Thereafter the Hamiltonian }i\ reads 



' Wo — 2 "-^i cos y (T3 



wi 



\/3/2 



(73) 



(74) 



At the second step we apply the transformation 



/3i /32 



-/52* 



(75) 
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where (3i, (32 are diagonal 2x2 matrices that must obey the unitarity condition 

|/3iP + |/?2p = l. (76) 



Supposing (3i^2 to be real and performing transformation (j75|) we come to the diag- 
onalization condition in the form 

M - Pi) + (Ai - X2)PiP2 = 0, (77) 

where Ai, A2 are 2x2 diagonal matrices and ^ is a parameter 

Ai = Wo — 3/2 (Ji cos 0(73, (78) 
cos 6 

A2 = -uJo-l/2uJi (T3 (79) 

cos a 

^ = u;i\/3/2sin6l (80) 

Assuming (32 = /?i where // is a diagonal 2x2 matrix as well we come to the 
following expressions for the matrix elements of /i 



fi, = ki + ^l + kl (81) 

where 

ki = -rrr^ = A^ — X2i. 
2? 

Finally we get for the matrix elements of Pi^2 

Pi = 1/2(1 + A;2)-V2 (^A;, + yTTfc^) ' (82) 
Pi = l/2|l + ^i^| (83) 

Now one can evaluate the connection 1-form. It is convenient to represent it as 
follows: 

M3/2 A^'-X 

A = iU^dU = Ad(l) = I #, (84) 

A/2/ 

where all matrix elements of A denote 2x2 matrix-valued blocks, U = U1U2U3 
and Ui are determined by (fTOl) . (173]) . (f75l) correspondingly. Here tilde denotes a 
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transposed matrix. After some algebra we get for the matrix elements of ()84p 



A'^ = -/Ji/32(3 - cosq)(T3 + - sina/?2(Ti/5i, (85) 

^3/2 = (a3/2 + fe3/2C^3 + C3/2Cri)#, (86) 

03/2 = ^ (3/3n - 3/?i2 + P21 COS a - P22 cos a) , (87) 

63/2 = ^(3/3fi + 3/??2+/3|^cosa + /3|2Cosa) , (88) 

C3/2 = sina/32i/522, (89) 

-4i/2 = {ai/2 + h/2<7s + Ci/2 0-i)d4>, (90) 

ai/2 = I (3/3|i - 3/3^2 + Pfi cos a - /3i2 cos a) , (91) 

61/2 = ^{Spl + Spi^+Plicosa + Pf^cosa) , (92) 

Ci/2 = -7: sina/3ii/?i2, (93) 



where d4> = ividt. Note that as A does not depend on time, the final solution does 
not require T-ordering. It should be also emphasized here that in the non-adiabatic 
case we discuss the term A3/2 contains non-diagonal terms that is not the case when 
the adiabaticity condition is held [58]. Now the solution of the problem takes a 
particular form of (j45p : 

^(t) = e-^'^"W e-^'^i*^'0(O), (94) 

Formula (|94|) solves the problem of the evolution control for the system under 
consideration. The resulting quantum gate is entirely determined by A and the 
evolution law of the magnetic field, i.e. by a contour in the parameter space. Of 
course it is always possible to choose the parameters so that A turns out to generate a 
2-qubit transformation that produces a superposition of basis states. For this reason 
the gate can be thought of as a universal one [T^. Of course, a suitable speed of 
the parameters evolution can not be reached by rotation of the sample as it took 
place in the experiment by authors of [7l]- Nevertheless it is clear that this manner 
of control is not principle and one could imagine a situation where the parameters 
evolution is provided by the controlling magnetic field by adding a non-stationary 
transverse component. It should be also noted here that general formulas ()85p do 
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not provide an apparent way to realize CNOT gate or another common 2-qubit gate. 
They just give the evolution law of the system provided that the external parameters 
vary as shown above. To knowledge of the authors other examples of computation 
of a conditional geometric phase that would be both non-Abelian and non-adiabatic 
are absent. To provide the gates of common interest one has to invent some special 
case of the parameters variation which makes the generic evolution operator more 
simple and transparent. This subject is out of the scope of this article. 

4 Conclusion 

The approach developed in [75 j- [77j is to be applied in the models where it is not 
possible to reduce the computation of the geometric phase to the case of 2- level sys- 
tem. Among those relevant to QC one can point out e.g. the model with anisotropic 
Heisenberg ferromagnetism where the exchange term in (jSip is determined by a ma- 
trix of constants Jab rather than by a single constant J. In this case the Zeeman 
terms H^'^ no longer commute with the exchange term Jab Sa Sb and to derive 
the expression for the geometric phase it is necessary to consider a more general 
case of 4-level system. The problem becomes more complicated also if the superfine 
electron- nucleus spin interaction must be taken into account. It is the case for the 
Kane model of silicon QC [79]. The effective dimension of the system's Hamiltonian 
is then 16. It is hopeless to attempt to obtain an exact analytic expression for the 
system's dynamics which should be investigated numerically (see e.g. |78] ) but it 
is nevertheless possible to derive an exact formula at least for the adiabatic phase. 

One more problem to be mentioned here is interaction with the environment. 
It can appear to be important not only for such issue as decoherence but it also 
can in principle contribute to the geometric phase. The simplest way to see it is to 
consider the model described by (jSOp - (1541) . If the external electromagnetic wave 
field can affect not only the qubits but the phonons as well. The phonon degrees of 
freedom can produce Heisenberg-Weyl-like geometric phase that can fill the sign of 
the spin projection due to electron-phonon term ()53p . (I54p . Some more complicated 
interaction between the qubits and the environment makes it necessary to compute 
the geometric phase for a system with the symmetry algebra which is larger than 
su(2) and cannot be reduced to the last one (in the sense of the phase derivation). 

Other field of application could be multi-beam optical schemes for quantum 
computations where several energy levels must be included in the scheme to provide 
two-qubit operations. Besides of some special cases [80] it can require more general 
methods for the geometric phase computation. 
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